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1 Introduction 



In recent years the physics of the standard model (SM) has moved to the investigation 
of flavor physics and in particular CP violation in the B meson decays. In the case 
of the neutral B decays, the mixing induced CP violation measurements involves a good 
knowledge of both the mass and width differences for which there is an extensive literature. 

Most of the analyses involving physics of neutral B mesons assume the validity of CPT. 
The B meson with their special properties provide unique opportunities of testing the 
CPT symmetry. Given the eventual need for an accurate description of the SM CP 
parameters, and/or tests for new physics, it is pertinent to also account for a possible 
violation of CPT invariance and its consequences. In this regard, it is noteworthy that the 
existing experimental limits on CPT violation are not very stringent and thus it becomes 
necessary to allow for this possibility in measurements of B decays. An argument for this 
was presented in P which questions the validity of CPT theorem for partons which are 
confined states. 

It is usually assumed that CPT is an exact symmetry which is hard to break. One of the 
first attempts to consider the (spontaneous) breaking of CPT came from string physics 
For instance, CPT may not be a good discrete symmetry in many higher dimensional 
theories of which the SM is an effective low energy description [H] . In addition, the existence 
of mixed non-commutative fields could also break CPT j3]. There are many tests for CPT 
symmetry and we refer to pi for a review on this subject. 

We present in this article a study of CPT violations that originate on the mass matrix 
and the analysis is model independent. We express the effects in terms of a complex 
parameter 6, whose presence modifies the expressions for width and mass differences and 
affects the time development of correlated B states. Precise measurements of the time 
development of the states are sensitive to the presence of CPT violation. However, such 
studies may require large number of events and hence it is preferable to have integrated 
rates. This motivates us to consider correlated decays of B mesons which are time ordered. 
The procedure presented here requires a time ordering of leptonic/hadronic events (which 
decay happens first) without demanding any detailed time development of states. One 
reason for following this time ordering is because it enables for a direct extraction of the 
width differences jH]. The same procedure is also expected to be sensitive to the presence 
of nonzero 6 values which can affect the width difference. We use the time ordered rates to 
form asymmetries which are sensitive to CPT violation. We find that our results depend 
on the width difference y = AT/ (2r) and thus for large CPT violating parameter ~ 0{y), 
the effects of CPT can become significant. Our analysis is more sensitive to Bg decays 
where the effects of CP violation can be neglected. We expect that the results obtained 
here will complement the existing tests for CPT violation, and in particular for the Bg 



2 



system where there is no information available as yet. For a recent general analysis of 
CPT violation in neutral mesons, and its extraction we refer to |2j. 

Briefly, the possibility of CPT violation has been studied in B decays using both time 
dependent and integrated methods. There are several analyses performed which in most 
cases also involve a detailed time dependent study and/or flavor tagging PIHI ITIHITH IT^ . 
In comparison with other known CPT studies, we find certain qualitative similarities as 
well as differences. Following our analysis, we find that we can restrict the strength of the 
CPT violating parameter to about 10% which is similar to an earlier detailed estimate jH]. 
We point out that if the strength of CPT violating parameter is below this limit (~ 10%), 
then using our approach, it is hard to pinpoint genuine CPT effects. In general, we find a 
sensitivity to CPT effects even if the width difference is vanishingly small. This we expect 
due to the particular time ordering procedure which can also extract corrections to the 
width difference due to nonzero 6. In general, in the presence of nonzero width difference, 
the CPT violating observables can measure both the real and imaginary parts of the CPT 
parameter; while on the other hand, for zero width difference, the effects are sensitive only 
to the real part of the CPT violating parameter ITU]. The latter feature is similar to what 
we find in our present analysis when we assume negligible CP effects. On the other hand, 
as observed in P^, for zero width difference, the CPT violating observable is sensitive to 
the imaginary part of the CPT violating parameter. This is different in our case (when 
we neglect CP violation), as we shall show when we discuss the Bg system. 

Our paper is organized as follows. In the next section, we present the general formalism for 
the time evolution of a neutral B meson in the presence of CPT violation. In section 01 we 
define the observables which can extract signals for CPT violation and are theoretically 
clean. This is followed by a brief numerical analysis for a specific decay channel and we 
compare the results with the B^ system. We finally conclude with a brief summary in 
section HI 

2 The time evolution 

2.1 Basic setup and definitions 

In this section, we review the basic formalism for the decay of correlated BB pair in 
the presence of CPT violation Starting with a quantum state which is a linear 

combination of the B and B states denoted as 

\m) = Mt)\B)+Mt)\B) , (1) 
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the time evolution of the BB system can be described by a two dimensional Schroedinger 
equation 



dt V Mt) J \M2i M22) V Mt) 
The mixing matrix is non-Hermitian and can be written as 

M = M-'-T, (3) 

with M and V being Hermitian 2x2 matrices (M = M^, V = F^). Invariance under CPT 
gives 

^11 = M22 Mn = M22 and Tn = , (4) 

while there is no constraint on the off-diagonal elements. As a result, to test CPT in- 
variance one needs to parametrize the difference of the diagonal elements of the mixing 
matrix. A sensible parameterization has to be invariant according to a re-phasing of the 
meson states ^1]. In particular, when a re-phasing of the type \B' ) e*^ \B ) is done, the 
anti-meson state is altered according to B'^ ^ e^*'^ )• Following this, we have in the 
matrix elements 

Mu ^ Mu ; M22 ^ M22 ; M12 ^ e-^'^'Mi2 ; -M21 ^ e^'^M2i . (5) 

This means the diagonal elements in Q remain the same after re-phasing while the product 
of the off-diagonal elements are re-phase invariant. Besides, the eigenvalues of Q are re- 
phase invariant. This leads to several possibilities to parameterize CPT violation and in 
this analysis, we choose a parameter 

In addition, there are other possible parameterizations which are rephase invariant and 
have been used for CPT studies. Among them we refer to parameterization in references 

[Biiiiin]. 

In the presence of 5 7^ 0, the width and mass differences of the two B states are obtained 
by calculating the eigenvalues of Q. We obtain these to be 



Ai = Mil + ^yMi2M2i 



^+4+2 



A, 



M22 - VMuM 



21 



52 



6 



(7) 
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We define A4ij = rriij — iTij/2 and 
Ai — A2 = 



2y/Mi2M2i\l— + 1 = Am - ^Ar . 



Therefore upon equating tlie real and imaginary parts in (jH| results in the width 
and mass difference 



Ar =ri-r2 = -4|mi2|ImWl- 
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— zRe 



mi2 V 4 



ri2 16^ 



— + 1 



Am = mi — m2 = 2 |mi2| Re^ / 1 — --^ — — 



\mi2\ 



mi2 V 4 



Furthermore, from the eigenvalue equation, 

MnPi,2 + Mi2qi,2 = Ai,2Pl,2 , 
M2lPl,2 + M22qi,2 = Ai,2gi,2 , 



we obtain the ratios 



Pi 

P2 



M 



21 



M 



12 



52 5 
^+4+2 



M 



21 



M 



12 



52 5 

1 + 



Introducing the mass eigenstates for the B mesons as 

1 

1 

the time development of B and B is evaluated to be 



\Bi) 

\B2) 



[pi\B) + q,\B)] 
[P2\B) + q2\B)] 



m)) 



l+UJ 



[U\B) + K,f^\B)] 



1 



m) = y--{u\B) + Kj^\B)) 



with 



/-t 
/- 



g-^^i* + cue" 



iA2t . 



= e 


-iXlt _j_ 


qi_ 


1 


pi 


~ ) 



-iX2t 



UJ 



P2qi , _ 
UJ = and UJ = UJ . 

q2Pi 



The information on CPT violation is encoded in the complex parameter 

uj = = X _|_ (5 _|_ — 1_ higher order terms in 5 . (15) 

q2Pi 2 

Given that 6 is re-phase invariant, u is also re-phase invariant and deviates from unity in 
case of CPT violation. 



2.2 The decay channel BB l^X^f{f) 

In the previous section we described the formalism for CPT violation that occurs in the 
mass matrix. We apply it to the decays of a correlated BB pair where one of them 
decays semi-leptonically while the other decays hadronically, i.e., BB l^X"^ f{f). The 
amplitude for a neutral B meson decaying into a final state with a lepton l^X^ at time to 
and the B decaying in to a final state / at time t can be expressed as 

A[/+(to), Rt)] = { X-t I B{t,) > ( / I B{t) > + C ( I B{t,) > ( / I B{t) ) , 

(16) 

with C denoting the charge conjugation of the BB pair. The individual decay amplitudes 
for the hadronic channels are 

{f\B) = A^e^'^'e'i' ; {f\B) = A2e''^^e'i^ , 

{f\B) = A2e-"^'e''^' ; { f \B) = Aie-^^'e'"' . (17) 

In (fTTjl . the Ai denote the absolute values for the amplitudes and rji and 0i are the strong 
and weak phases respectively. Using this, along with the notations 

A p2i(t>m 



Ai + Re[6] 

with the phase due to mixing denoted by the amplitudes are now: 
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I i + u; I 



[Mto)f.it) - /-(to)/+(t)] + ^e^^e^'' [/+(to)/+(t) - f-ito)f-it)] 



[f-{to)f-{t) - U{to)U{t)] + ^e^*e^''[/-(to)/+(t) - /+(to)/-W] 

{ [f-ito)Uit) - Uit,)f4t)] + r lel e-*e^^ [/-(^o)/-(t) - fMUit)] } ■ 

(19) 

In fjTTH) . F denotes the amplitude for the decay i? ^ l^X^ and the corresponding amphtude 
for the CP conjugated process is denoted by F. The decay to a lepton is measured at to 
while t is the time when the decay to state / occurred. An explicitly time dependence 
in the l.h.s. of ()19p is assumed. The BB pair is taken to be the one produced at a T 
resonance and hence, the BB charge parity, C = — 1 in In addition, the state / is 

chosen to be a CP odd eigenstate. This sets 77 = and r = —1. 



3 The CPT violating rate 

Following the above construction, for the decay process BB — ^ l^X^f we define two time 
correlated observables: (i) Rs denotes the number of events in which the hadronic decay 
precedes the semi-leptonic one (which in our context the time ordering is t < to) and 
(ii) similarly, we define the number of events where the semi-leptonic decay precedes the 
hadronic decay denoted by Rl. To illustrate, if we choose positively charged leptons, we 
have the rates 

00 to 00 00 

Rs = Jdtoj dt\A[l+f]\' ; Rl [t f] = j dto j dt\A[tf]\'' . (20) 

00 io 
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Using (j2ni), for CP eigenstates ^ we can define the following ratios 

Rs [l+f]±RL[l-f] 



R^ 
Ri 



R[l+ f] + R[l- f] 

rl [i^n-Ri [i-f] 

R[l+f] + R[l-f] 



_ Rs [/+/] - Rs [l-f] .... 
- R[l^f] + R[l-f] ■ ^21) 

where the rates without the subscripts L or 5 denote total time integrated rates without any 
time ordering. The exact expressions for the above ratios can be obtained by substituting 
the expressions for the rates given in the appendix (sectional) of this paper. Due to the 
extensive nature of the analytic expressions involved, we avoid from presenting them in 
the main text as they are also not illustrative. 



3.1 Numerical analysis 

The effects discussed in this section were calculated using the exact expressions for the event 
rates which are given in the appendix (section^ of this article. In these expressions and for 
our specific case, where we are interested in CP odd eigenstates, we set rj = and r = — 1. 
In this analysis, we have varied both Re[5] and lm[S] in the range between —0.5 to 0.5. It 
has been shown in ^2] that this range is not excluded yet by the of the recent Belle data, 
according to which \mB — rn§\ ~ lO^^^m^ J7]. Figs C] and |2 show the results for the Bg 
system and for appropriate values of the variables. We mainly analyse the Bg system and 
also briefly present the results for the Bd system where one has to account for CP effects 
arising from mixing. 

In Fig. ^we show the ratio R~[ as a function of Re [5] and for three values of y consistent 
with the current estimates [TH]. In the CPT conserving limit, the value for R^ can be read 
off on the axis Re[S] = and its range is —0.10 < Ri < —0.025. When Re[6] ^ there 
are deviations from this range and for Re [5] ~ ±0.1 the deviations for R^ is ~ 0(0.10) 
from the central line. Thus if a deviation of that order is observed it will come either from 
violation of the CPT symmetry or there is a large width difference of y > 0.12 either of 
which is very interesting. 

Taking y = 0.12 and x = 19, in Fig. |21we show the ratio i?2 as functions of Re[5] and 
Im[(5]. Again there are sizable deviations when Re[5] ^ 0. The structure of the curves in 
Fig. ^and Fig. [21 can be accounted for examining the small 6 limit of the ratios in (pT|. 

^This choice sets f — f in H16|) . 
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For instance, defining y = Ar/(2r), x = Am/T, in the small 6 limit which we take such 
that (Re[5],Im[(5]) < 0.1, we find 



R- ^ | + Re[5] 



2 (a;2 + 11 



+ 0(6^). (22) 



Thus, a linear dependence of Ri as a function of Re [5] is evident near the origin of the 
plot in Fig. H Note that in the limit y = and in the absence of any CP or CPT violation 
we would expect an equal number of and /~ events. 

Similarly, we observe in Fig. |21that for small Re [5] a linear behavior in R2 and this is again 
understood by examining the R2 in the small 6 limit wherein 

R2 ^ ^(T^(f^^[^]^ + M5]^(2/ + 1)) ^^fory«x. (23) 

In Fig. 121 for all values of Re[5], the sensitivity to lm.[6] is almost negligible. This feature 
is shown by the fiatness of R2 as a function of Im[5]. In addition without imposing the 
small Re [5] limit, we have also numerically checked the dependence of on lm[6] and we 
find it to be too small ~ 0(0.001). This can be attributed to the largeness of x for the Bg 
system which is also brought out in our approximate analytic expressions discussed above. 

As a consistency check, upon using (j22I) and we easily find the relation 

i?r - 2/?2 = I . (24) 

This relation can be verified by combining the results in Fig. ^ and Fig. El where we 
find our numerical results are consistent with (j^^ for the small 6 limit. Furthermore, as 
is evident, the ratio R2 is nonzero only in the presence of CPT violation and serves as 
consistency check for (j^ . In addition, we note that in the case of system, for y <ti x, 
in the small S limit, we have the approximation 

R - R - + ~ (25) 

R2-Rz-^ + — ■ (25) 

This calls for the following remark. As a consequence of the equality, the sum of these two 
ratios {R2 + -R3) would prove to be a good observable as it does not require time ordering. 
However, since the equality is a result of a specific limit, we may still plot them separately. 

It is interesting to ask what is required in experiments in order to observe these CPT vio- 
lating effects. First of all, the effects which are discussed here are for C = — 1 state, which 
is satisfied in e^e~ colliders; provided these machines are tuned at higher energies to be 
able to produce B^Bg pairs. A promising channel is B^B^ — > l^X^J/ipcj) where the lepton 
and J/ip'P are detected at two different times with the identification of the time ordering 
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at the decays, i.e. which decay happened first. The measurement of the ratios requires an 
accuracy of 1% and we need about 10^ decays to the channels under consideration. Given 
the semi leptonic branching fraction of Bs to be 10% and the decay to J /ipcf) to be 10^^ we 
have a branching fraction for this process of ~ 10~^. Consequently, we require roughly 10^ 
BgBs pairs produced through a C = — 1 state. This also suggests that with the comparable 
number of events, a reasonable bound can be set for Re[5]. Alternatively, if no significant 
deviation from the above mentioned range for say is found, then one can conclude that 
the data are consistent with zero or rather small CPT violating effects. Here, by small, 
we mean Re [5] < 0.1. This limit can be improved if in future we have better handle on the 
values of y which currently suffers from large errors. 

We now turn to discuss briefly the effects of CPT violation for the Bd system and we choose 
the process, BdBd — ^ l^X^ J /ipKs- In contrast to the Bg mesons, for the case of the Bd 
system, one cannot neglect the presence of the CP violating phase. We show in Figs. El 
and El the variation of Rf with Re[(5] and $. In these two plots, we have set Im[(5] = 
while setting x = 0.755 and y = 10~^. The variation of R]^ with Re[5] is significant and 
we show this in Fig. El The width in the plot is due to the allowed range of values in x 
and $. Presently, their ranges are sin<l> = 0.735 ± 0.055 [W with x = 0.755 ± 0.015 and 
y ~ 10~^ The range of is now smaller as compared to the Bg case and will require 
more events to reach the bound of Re[5] < ibO.l Finally, the ratio R^ depends on both 
Re[5] and lm[6] and illustrated in Fig. El For this plot, we use the values sin$ = 0.735, 
X = 0.755 and y = 10~^ . It is interesting that, we can observe the impact of lm[6] which 
was not so significant in other plots. In contrast to this analysis, the effects of lm[6] (with 
Re [5] = 0) on width measurements and CP asymmetries for the B^ mesons has been 
examined in reference |12j . 

Finally, we address the question of nonzero CP violation. Although, the present analysis 
is most effective for the Bg system, we note that the ratios in (f?T|) can still be used to 
test for CPT violating effects even if the system has a nonzero CP violation. In the pure 
CP conserving limit, for instance, independent of y, the ratio Rf can show a numerical 
deviation from the value of 1/2 in the presence of CPT violation. This becomes evident 
by examining the small S limit where 



In the presence of CP violation the ratio Rf gets modified as follows, Rf —>■ Ri+r'^, where 
the typical strength of r+ is a linear combination of the form r+ ~ (sin$/a;, cos $ -Re [5] /x). 
Thus, for the Bg system, given the expectations for x ~ 19 and a small CP phase such that 
sin$ ~ 10~^ the ratio Rf is sensitive to CPT effects for all values of (Re[5], Im[(5]) > 




^ + Re[5](|)-Im[<5] 



2(1 + ■ 



(26) 
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In the case of CP and CPT conservation we expect the rates Rsll^^f] = -Rl[/^/] = 
1/2{R[1~ f] + R[l~^ f]) which means that the decays into hadrons before the leptonic de- 
cay are equal to the hadronic decays which occur after the leptonic decay. The equality 
of the above decay rates is modified in the case that there is CP or CPT violation. For 
Bfi decays the CP violation has been observed and will produce such a difference; an in- 
teresting question here is whether the CP-violation present in Rf is consistent with the 
CP violation established in other channels, like J/ipRg. Deviations will be attributed to 
new physics and/or CPT violation. Furthermore, in light of the recent discrepencies in 
the measurements of CP violation in the decays —>■ J/ipKg and —>■ (pKg one could 
envisage new physcis contribution involving new penguin operators j22j. Such operators 
could induce modifications to mixing matrix defined in and hence modify the expected 
signals from the ratio R^. In the present analysis, we do not consider these possible cor- 
rections which may be required if the present signals for new physics persists and becomes 
statistically significant. The effects in B^ decays will be even more interesting, because 
it is still possible that in extended models sin($)/x is larger than 0.04 and comparison 
among various channels, like J / tpKg with (pKg will be interesting to establish the origin of 
the effect; whether the new physics originates in the CP or CPT sector or from both. 

However, we note that in principle, by strictly measuring Rf alone, one cannot ensure 
the deviation from 1/2 as an unambiguous signal for pure CPT violation or pure CP 
violation. On the other hand, any corrections due to CP can affect the ratio i?f only 
at the sub-leading level which is suppressed by (Re[(5], Im[5]) and hence is not significant. 
In the present numerical analysis, we consider much larger values (Re[(5], Im[5]) ~ 0(0. 1) 
which are not excluded by current available data. Clearly, at this level the effects due to 
CP violation are expected to be small for the B^ system. In the case for the B^ system, 
the impact of $ 7^ can be important and must be studied explicitly in the analysis of the 
data. 

4 Summary 

We discussed indirect CPT violation as it appears in the B meson system. The break 
down of the CPT symmetry occurs in several theories and modifies the time development 
of the states. The development is characterized in terms of a parameter 6 which is phase- 
convention independent Consequences of the break down of the CPT symmetries 
manifests itself in the time development of the states. Our attention is addressed to the 
production of BgBs and B^Bd pairs in a odd charge conjugation state. In the case of B^ 
and Bg system, the decay depends on the CPT violating parameter 6 while the influence 
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of the CP phase is neghgible. One way to observe the difference (due to 5) is to study 
the time development of the decays. However, because of the large number of the events 
required, we propose time integrated rates. Once the BgBs pair is produced, we encounter 
two different decays; one of them can be semi-leptonic and the other one can be hadronic. 
We also defined a time-ordering among them, meaning which one of these events occurs 
first. 

The time ordered ratio Rs denotes all the events where the semi-leptonic decay follows 
the hadronic one; similarly, Rl includes the events where the semi-leptonic decay occurs 
before the hadronic jO] . These two rates are different for two reasons: the width differences 
and because of CP and/or CPT violation. The ratios were calculated using the formu- 
las given in the appendix and the results were presented in section 13.11 We calculated 
the asymmetries involving Rl,s as functions of the parameter 5 and the observed width 
difference y. We summarised our numerical results in section 13.11 and through Figs. Q to 
ini The formalism presented here holds for Bd system as well and our conclusion is that 
experiments with 10^ B^Bg or B^B^ pairs will be able to restrict the Re [5] smaller than 
10% or otherwise observe an effect. These considerations are within reach of the present 
experiments [THES]- 

It is important to note that due to the specific nature of the time ordering, the various 
ratios considered here are sensitive to y and hence simultaneously to any CPT pollution 
which can affect y. The impact of CPT violation was shown qualitatively for the small S 
limit. In this limit, for zero CP violation, which is a good approximation for the Bg system, 
the ratios show a very simple dependence on the parameters x and y; thus allowing for a 
direct dependence on 6 under suitable limits {x ^ y). We remark that for this limit, the 
impact of CPT effects becomes sensitive to the strength of Re[5] and Im[5], more so to 
Re[6] due to large x ^ y; and thus tests the magnitude and phase of the CPT violating 
parameters. To get a qualitative feeling for the CPT effects, we find from ()22|) that for 
Re[5] ~ y, one can envisage the ratios to exhibit CPT violating effects ~ Re[5] and thus 
can be large^ ~ 0(0.10). In a similar spirit, from we can expect deviations (here of 
0(0.02)) which is relatively independent of the parameters x and y. To our knowledge, 
this situation is in contrast with many other interesting alternative methods known in the 
literature where such a simple dependence is perhaps not observed in the limit of a small 
CPT violating parameter and thus show a different sensitive to CPT effects than the 
method prescribed here. 

^As mentioned earlier, this effect is observable provided Re[(5] does not get washed out due to the errors 
in y. 
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A Appendix 



In this section, we present exact expressions for the rates calculated from (|Tn|l using the 
definition 



oo to 



RsPifJ)] = jdtoj dt\A[l^{fJ)]\' , 





oo oo 



RiPifJl] = jdtoj dt\A[l^{f,f)]\' . (27) 



to 

In our notation, Q represents the phase of uj ^ 1 + 6. Defining the overall prefactor 



^ 2r2|cj|2(l+x2)(y2_l) ' (^^^ 

we have: 

Rs[l+f] = -ir|2(x2 + y2)|^|2|^|2 + 2|ea;|r[(l + |/)(x2 + y)|cu|cos(r7 + $) 

- (y - l){y - x^) cos(?7 + $ - fi) + x(l + y)\\u)\ sin(r7 + $) + sin(?7 + $ - f])] 
+ r2[(l + x2)(l + |c^|2 + y(|c^|2- 1)) + 2(2/2- l)|o^|(xsinl]-cosl])] 

i?L[/+/] = ir|-2(x2 + y2)|er|c.p + r2(l+x2)(-l-y+(-l+y)|c.p) 

+ 2r|cj||(-l + y^){rco^Vt - x\i\{\uj\ + cosf2) sin(?7 + $)} 
+ (1 + 2/){rx(-l + ?/) + [x^ + y)\i\ •mviji + $)}sinfi 

+ |^|cos(?7 + $){(l+y)(x2 + y)cosfi + (?/-l)[(x2-y)|cj|+x(l + ?/)sinn]}|| . 

(29) 
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Rs [l-f] = X| - 2r2(a;2 + y^) + (1 + - y + {y - l)\io\') - \^\ [2r ((1 + y){x^ + y) 

X cos(77 + $) + (?/- - y)\uj\ cos{r] + ^ - Q) + x{l + j/)(sm(r/ + $) 

+ |a;|sm(?7 + $-fi))])] + 2|^|2 (y^ - l) |a;| [cosfi + xsinQ] | , 

Rl [l-f] = -Ki^2r'[x' + y'] + (1 + x^jl^ril + + |/(|c.|2 - 1)) - 2{y' - mm^\ ^osQ 

+ rx{l + \uj\ cosVl) sm{r] + $) - 2(1 + y)\^Lu\{-x{y - 1){ + r(a;^ + |/) sin(?7 + $)) 

X sinQ + 2r|^| 003(77 + '^)(-(l + Z/)(2^^ + v) 



X |a;| cosf2 + (y — l){—x'^ + y + x{l + yjlto'l sinf2)) 



(30) 



Rs = -K^2r'(x' + + (1 + x^)(l + + y(kr - 1)) + 2(y' - l)|a;| 

X (x sin Q - cos n) + \u;\ [2r\^\ ((1 + y)(x^ + y) |<^| 003(77 - ^) + (v - - v) 

X cos(77 - $ + n) + a;(l + y)(|a;| sin(77 - $) + sin(77 - $ + O))))] > , 



Rl [l^f] = - (1 + a;2)(l + y) + ((1 + a;^)(2/ - 1) - 2r\x^ + 2/^)|er)|a;p + 2r(a;2 - y) 
X (y - l)|e||c^|'cos(77 - $) + 2(y2 _ cosQ + 2\u\ [r(l + y)(x2 + y)\^\ 
X cos(77 - $ + Q) + x(y^ - l)(sinQ + r^(|a;| sm(77 - + sm(77 - $ + Q)))] I . 

(31) 
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Rs [l-f] = X| - 2y' - r^ier + r'\^\\y + - y\u^\') - x'{2 + r'\^\\l + y + \cu\' - y\u;\')) 

- 2r(l + y){x' + y^)\^\ 005(77 - $) + 2r{y - l)|e|[r(l + y)\^u\ 

X cosf2 + {y — a;^)|a;| cos(?7 — $ + x{l + y){sm{rj — $) + r|^|a;sinf2) 

+ 2rx{y'^ -l)\^uj\sm{r)-^ + n)> , 



RL[l-f] = -X^ 2(^2 + ^1^1 ^|^|((i + ^2^(i + |^|2 + ^(|^|2_i)) + 2|a;|(y2_i) 



X (-cosQ + xsinQ)) + 2 (7/ - l)(y - x^) 003(77 -$)- (1 + + y) 



X |a;| cos(77 - $ + Q) + x(y^ - l)(sin(77 - ^) + \uj\{r] - ^ + fi)) 



(32) 
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Figure 1: Ratio Ri versus Re [5] for 
the Bg system. The three curves cor- 
respond to \y\ = 0.05 (upper line), 
\y\ = 0.12 (middle line) and \y\ = 0.19 
(lower line). 



0.4 



Figure 2: 3D Plot for R2 showing the dependence on Re [5] 
and lm[6] for the Bg system. 
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Figure 3: The ratio Rt as a function of $ t^. a r^^ mr ^■ 

^ figure 4: i iie same as ig. |o| tor the ratio 

for the Bd system. The phase $ is given in 



radians. 




Figure 5: Dependence of on Re [5] for 
the Bd system. The width in the plot is 
due to the errors in y,x and $. The range 
of these values are specified in the text. 




Figure 6: 3D Plot for the ratio R3 its de- 
pendence on Re [5] and lm[6] for B^ system. 
One can see that the imaginary part is de- 
tectable. 
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